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Abstract. We construct and discuss Hadamard states for both scalar and Dirac spinor fields in a large 
class of spatially flat Friedmann-Robertson- Walker spacetimes characterised by an initial phase either 
of exponential or of power-law expansion. The states we obtain can be interpreted as being in thermal 
equilibrium at the time when the scale factor a has a specific value a — a^. In the case ciq — 0, these 
states fulfil a strict KMS condition on the boundary of the spacetime, which is either a cosmological 
horizon, or a Big Bang hypersurface. Furthermore, in the conformally invariant case, they are conformal 
KMS states on the full spacetime. However, they provide a natural notion of an approximate KMS state 
also in the remaining cases, especially for massive fields. On the technical side, our results are based 
on a bulk-to-boundary reconstruction technique already successfully applied in the scalar case and here 
proven to be suitable also for spinor fields. The potential applications of the states we find range over 
a broad spectrum, but they appear to be suited to discuss in particular thermal phenomena such as the 
cosmic neutrino background or the quantum state of dark matter. 

Contents 

1 Introduction 2 

2 Prom the Geometry to the Classical Field Theory 4 

2.1 Friedmann-Robertson- Walker Spacetimes with a Lightlike Cosmological Boundary 5 

2.2 Classical Scalar Fields on FRW Spacetimes 8 

2.3 Classical Dirac Fields on FRW Spacetimes 10 



1 



3 Quantum Field Theory in the Bulk and on the Cosmological Boundary 16 

3.1 The Bulk-to-Boundary Correspondence for Real Scalar Fields 16 

3.2 The Bulk-to-Boundary Correspondence for Dirac Fields 18 

4 Boundary States and their Hadamard Bulk Counterpart 23 

4.1 The Scalar Field Case 24 

4.2 The Dirac Field Case 27 

4.3 The Bulk Remnant of the Boundary KMS Condition 29 

4.4 Relaxing the KMS Condition on the Boundary 31 

4.5 Approximated Thermodynamics 32 

5 Conclusions 33 
A On the Definition of Dirac Fields 34 
B Perturbative Construction and Analysis of the Classical Solutions 36 



1 Introduction 

The mathematically rigorous formulation of quantum field theory on globally hyperbolic curved 
spacetimes has undoubtedly witnessed terrific progress during the past years. Our understanding 
of the quantization scheme within the algebraic framework initially allows us to specify which 
quantum states are physically relevant for any free field theory; based on these, one can then 
construct both Wick polynomials and interacting field theories on non-trivial backgrounds in a 
perturbative manner [BrFrOO, HoWaOl, HoWa02, BFV03, HoWa03, BrFr09]. 

Hence, one of the main challenges is to apply this mathematical scheme to concrete problems, 
and the first and foremost examples are certainly found in the realm of Cosmology. Accord- 
ing to the standard modern paradigm, the Universe can be described by a homogeneous and 
isotropic solution of Einstein's equations, a Friedmann-Robertson- Walker (FRW) spacetime. 
This requirement entails that the metric is fully determined by a single dynamical quantity, 
the so-called scale factor a(t), and a parameter which fixes the topology of the constant-time 
hypersurfaces to be spherical, flat, or hyperbolic. In view of recent measurements of the evo- 
lution of the Hubble function, we only consider the case of flat spatial sections - so-called flat 
FRW spacetimes - in the present work. Furthermore, we demand that the scale factor of the 
FRW metric displays an inflationary behaviour of either exponential or power-law type at early 
times. Our aim is to show that, under these circumstances, it is always possible to single out a 
so-called Hadamard state for both scalar and spinor quantum fields. The Hadamard condition 
is of particular relevance since it guarantees that the UV behaviour of the state mimics the one 
of the Minkowski vacuum, which in turn assures that the quantum fluctuations of observables 
such as the smeared components of the stress-energy tensor are bounded. 

The scenario of a free scalar field has been extensively studied in the past and it has been 
found that, on a FRW spacetime with flat spatial sections, a Hadamard state invariant under 
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the metric isometries can be singled out in various ways, see [JuSc02, O107, Kiis08, SchVe08] and 
[DMP07, DMP09]. Here, we will focus on the approach first introduced in the last two last cited 
papers, where the notion of a distinguished asymptotic Hadamard ground state was discussed by 
means of a so-called bulk-to-boundary approach. This is a procedure which calls for identifying 
a preferred null, differ entiable, codimension 1 submanifold of the full spacetime - the boundary 
- on which it is possible to construct a *-algebra of observables which contains the bulk one 
via an injective *-homomorphism. The merit of this construction is that the peculiar geometric 
structure of the auxiliary submanifold allows for identifying a boundary ground state whose 
pull-back in the bulk is both invariant under isometries and, more importantly, of Hadamard 
form. 

The first aim of this work is to extend this procedure in order to obtain additional physically 
interesting states for free fields. As first realised in [DMP10], the structure of the boundary is 
such that it is also possible to naturally construct states for the boundary algebra of observables 
which are thermal, that is, they fulfil a KMS condition with respect to a suitable boundary 
translation. We will show that this idea can be implemented successfully in a cosmological 
framework and for scalar fields, yielding states which are both Hadamard in all cases and con- 
formal KMS states in the case of massless fields conformally coupled to the scalar curvature. 
More interestingly, it turns out that, even though the KMS condition is not exactly fulfilled 
for massive theories, approximate thermodynamic relations always hold at early times. Hence, 
the states we find can be interpreted as being asymptotic thermal states in the early universe, 
where the departure from strict thermal equilibrium is controlled by the magnitude of the scale 
factor times the mass. We generalise these results even further by introducing Hadamard states 
which fail to satisfy a KMS condition on the boundary, i.e. at a(t) = 0, but rather fulfil an 
approximate KMS condition at a(t) = ao > 0. 

The second main goal of this paper is to prove that the bulk-to-boundary procedure can be 
extended to the case of free massive Dirac fields. We show that it is possible to construct a 
suitable boundary algebra of observables encompassing the bulk one also in this case, and that 
one can naturally identify preferred states on this boundary algebra whose pullback is always of 
Hadamard type and either a strict conformal KMS state for massless fields, or an approximate 
KMS state in the case of non-vanishing mass. 

We want to stress that our analysis is not performed only out of a mere mathematical interest, 
but that there actually cogent physical reasons to study the topics we discuss in this work. On 
the one hand, one of the most interesting applications of the algebraic quantisation scheme 
to Cosmology is the analysis of solutions to the semiclassical Einstein equations for massive 
fields. Indeed, it has been found in [DFP08] that a scale factor driven by a massive scalar field 
in a Hadamard state naturally evolves to a late time de Sitter solution, thereby offering an 
explanation for dynamical dark energy. It is certainly mandatory to analyse this phenomenon 
also for fields of higher spin and to see if more observed features of the cosmological evolution 
can be modelled by the sole effects of quantum fields in more general Hadamard states. In fact, 
the very states constructed in this paper have already been analysed in this manner to confirm 
that the backreaction of free Dirac quantum fields on the spacetime curvature is qualitatively of 
the same nature as the one of free scalar fields, and it has been found that the energy density of 
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dark matter, usually considered as an ad- hoc classical quantity, can be accurately modelled by 
free quantum fields in an approximate KMS state of the kind introduced in this work [DHMP10]. 

On the other hand, an important object of study in modern cosmology is the cosmic neutrino 
background, which is usually considered to have a temperature of approximately 1.9 Kelvin (see 
for example [KoTu90]). However, since it is now well-known that neutrinos are massive and, 
hence, do not fulfil a conformally invariant equation of motion, it is impossible to associate to 
them an exact conformal KMS state. Therefore, the approximate KMS states for massive Dirac 
fields we introduce seem to be natural candidates to discuss the cosmic neutrino background 
accurately in terms of quantum field theory in curved spacetimes. 

The paper is organised as follows: section 2 starts with a recollection of the geometrical 
properties of the spacetimes we consider and closes with a study of the properties of both a 
classical scalar and a classical Dirac field theory living thereon. It is particularly emphasised 
how to solve the Dirac equation in a flat Friedmann-Robertson- Walker spacetime in terms of a 
vector- valued, diagonal, second order hyperbolic partial differential operator. Section 3 instead 
is focused on the construction of the bulk-to-boundary correspondence. After recollecting how 
this mechanism works for scalar fields, we continue to fully develop it for spinors. Particularly, 
we prove the existence of an injective *-homomorphism from the bulk field algebra into the 
boundary counterpart. Section 4 is the core of the paper and it is here shown how to assign 
distinguished states for a free field theory on the boundary which are satisfying an exact KMS 
condition. Afterwards, the bulk counterpart is constructed, and it is shown that both in the 
scalar and in the Fermionic case the outcome is a Hadamard state which is in addition ther- 
mal when the classical dynamics of the bulk field is conformally invariant. In case conformal 
invariance is broken, we show that the result still fulfils approximate thermodynamic relations 
and that this even holds for states which are not exact KMS states on the boundary, but fulfil 
an approximate conformal KMS condition at finite times. Section 5 contains the conclusions, 
while the appendices are focused on reviewing the definition of Dirac fields on a four-dimensional 
globally hyperbolic spacetime and on proving the fall-off behaviour at infinity of classical Dirac 
modes in the considered flat FRW spacetimes. 

2 From the Geometry to the Classical Field Theory 

This section is bipartite: in the first part, encompassing the next two subsections, we recollect 
some facts already discussed in [DMP07, DMP09, PinlO]. Particularly, we define the class 
of spacetimes we are interested in and recollect their geometric properties. Furthermore, we 
sketch the behaviour of a classical real scalar field living thereon. In the second part, we focus 
on classical Dirac spinors and provide a detailed account of their dynamics on cosmological 
spacetimes. 
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2.1 Friedmann- Robertson- Walker Spacetimes with a Lightlike Cosmological 
Boundary 

In this paper, we shall consider spacetimes M as being four-dimensional Hausdorff, connected, 
smooth manifolds endowed with a Lorentzian smooth metric g whose signature is (—,+,+,+). 
Particularly, we shall be interested in homogeneous and isotropic solutions of Einstein's equations 
with flat spatial sections. Thus, the metric qfrw has the so-called Friedmann-Robertson- Walker 
form: 



which is here written in spherical coordinates. The coordinate t, also known as cosmolog- 
ical time, ranges a priori over an open interval ICR which we later constrain, whereas 
a(t) G C°°(I,M. + ), R + being the strictly positive real numbers. If one introduces the so-called 
conformal time r out of the defining differential relation dr = a~ 1 (t)dt, the metric becomes: 



which is manifestly a conformal rescaling of the Minkowski one. 

This rather simple consideration initially prompted our attention towards this class of space- 
times, since the flat background is the prototype of a large class of solutions of Einstein's 
equations with vanishing cosmological constant, namely, the asymptotically flat spacetimes (see 
[DMP05] or [Wa84] for a recollection of the definition). From a geometrical point of view, these 
are rather distinguished since, by means of a conformal compactification process, they can be 
endowed with a notion of (conformal) past or future null boundary, usually indicated as 9 ± . 
This rather special structure was successfully used in [DMP05] (see also references therein) as a 
tool to set up a bulk-to-boundary correspondence ultimately yielding the possibility to construct 
a quasi-free invariant Hadamard state for a massless scalar field conformally coupled to scalar 
curvature. 

Therefore, it seemed natural to wonder if a similar bulk-to-boundary correspondence and 
related results could be obtained also for FRW spacetimes. To this avail, the first question 
to answer is, under which conditions on the scale factor a(t) past or future null-infinity can 
be genuinely associated to (M, gFRw)- Indeed, if we require that the functional form of a(t) 
and its domain I 3 t are such that the domain of a(r) includes (— oo, — tq) ((to,oo)) for some 
ro < oo, then one can meaningfully assign to (M,gpRw) conformal past (future) null infinity 
by means of the conformal equivalence with Minkowski spacetime. As we would somehow like 
to consider the field theoretical constructions on null infinity as "initial conditions" rather than 
"final conditions" for the field theories in the bulk spacetimes, we shall restrict attention to 
cases where (M, gFRw) possesses a meaningful notion of past null infinity. However, all our 
constructions and results can be trivially extended to FRW spacetimes with future null infinity. 
Four simple examples of cosmological spacetimes with past null infinity are: 



ds 2 = -dt 2 + a 2 (t) [dr 2 + r 2 {dd 2 + sin 2 6dip 2 ] , 



(1) 



ds 2 = a 2 (r) [-dr 2 + dr 2 + r 2 (d9 2 + sin 2 9d<p 2 )] 



(2) 



n 



n 



a 



) a(t) = (^) 



n > 0, 7 > 0, t G I = (-oo,0) =>■ a 



( T ) ~ (-r(n+l)) 



n+l 
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b) a(t) = i, 7 > 0, t e I = (0, oo) =► o(t) = 



c) a(t) = (i)", n > 1, 7 > 0, t € / = (0, oo) a(r) = (^i)) 

_^ 

d) a(i) = eT,7>0, t€ I = (-oo, oo) a(r) = ^ 

Since only the behaviour of a(t) at the lower bound of / 3 t is essential, one can immediately 
obtain a large class of flat FRW spacetimes which posses past null infinity by requiring that 
the functional behaviour of a(t) is asymptotically of one of the above mentioned kinds in the 
early past. The FRW spacetimes one obtains in this manner from the cases b), c), and d) are of 
special cosmological relevance, as they describe a cosmic history characterised by an early stage 
of either de Sitter inflation - d) - or power-law inflation - b) and c). This is the physical 
motivation why we shall restrict to these cases in the following. The technical motivation is the 
fact that, as found in [DMP07, PinlO], in these spacetimes it is possible to establish a bulk-to- 
boundary relation for massive field theories, while, in FRW spacetimes which are asymptotically 
of class a), this is only possible for massless theories. To understand the reason behind this in 
simple terms, we briefly describe how past null infinity can be attached to (M, gFRw) if the 
above described requirements are met. 

Switching to the coordinates U = tan _1 (r + r) and V = tan~ 1 (r — r), the FRW metric reads 



ds 2 



a 2 (r(U,V)) 
cos 2 U cos 2 V 



. dUdV + Sin( ^ V \ d0 2 + sin 2 e V) 



(3) 



where the factor in the square brackets is the line element proper to the metric of the Einstein 
static universe (Me^e)- As a result, flat FRW spacetimes can be conformally embedded into 
this spacetime, and the locus U = — ^ in Me - corresponding to the "boundary" r+r = — oo of a 
FRW spacetime with the domain of a(r) containing (— oo, —To) - can be meaningfully considered 
as a smooth hypersurface 9~ diffeomorphic to Ex S 2 , i.e. past null infinity. A bulk-to-boundary 
construction as the one we seek to employ always requires to extend the definition of a field in 
the original spacetime - here {M,gpRw) - to a field on the boundary, namely, How this 
can be achieved crucially depends on the behaviour of Qfrw on past null infinity. In spacetimes 
which are asymptotically de Sitter, one finds that the prefactor of the square brackets in (3) 
becomes constant on Hence, in this case, (M,gpRw) can be extended beyond and the 
bad behaviour of gFRW on Q~ turns out to be rooted in an unlucky choice of coordinates. In this 
case, turns out to be a horizon of the FRW spacetime. In contrast, in case a), the prefactor 
in (3) diverges on while it vanishes in the remaining cases, i.e. in spacetimes with an early 
power-law inflation. Accordingly, 9~ describes the Big Bang hypersurface in the power-law 
case. On the field theoretic side, the failure of a l^'^L to be finite on is inherited by 

' cos^ U cos z V J 

solutions of the classical field equations and forces us to introduce suitable conformal rescalings 
of the field in order to be able to map a theory in a finite manner from (M, gFRw) to As a 
result, a field with mass m becomes a field with mass 



a 2 (r(U,V)) 

277 2~T7 m - 

cos^ U cos^ V 
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Its is now clearly visible why a bulk-to-boundary construction for massive fields can only be 
established in spacetimes which are asymptotically of de Sitter or power-law inflationary type. 

For the technical details of our analysis it will be convenient to discuss all constructions on 
the level of conformal time r. Hence, let us recall the related requirements on the scale factor 
a(r) of the flat FRW spacetimes we wish to work with in the following, viz., 

l l+5 1 A da{r) (l + 6h 1+S 1 \ f ^ n . AS 

<r) = ^ + 0[j-^y ^ = _i_^ + O ^-^j,*>0. (4) 

Here, the formulae are meant to hold asymptotically for r — > — oo. Furthermore, e is a small 
positive constant that has to be added in order to obtain the desired regularity for the classical 
solutions near null infinity in the case 5 = 0. In favour of notational simplicity, we consider the 
special case a(t) ~ t^ 1 <^ a(r) ~ e 7 " 7 to be represented by the limiting case S — > oo in the 
above formulae. 

As the interplay between bulk symmetries and the vector d T with the geometric structure 
of S~ will play a pivotal role in the following, we shall now recall a few more details on the 
geometry of 9f — . In the asymptotic power-law case, 3 s- has manifestly the same properties as the 
past conformal boundary of an asymptotically flat spacetime, extensively discussed in [DMP05] 
and the references therein. In the asymptotic de Sitter case, the structure of 9 _ is different due 
to the fact that it is a horizon, a full discussion can be found in [DMP07]. In more detail, setting 

9 

= 2 cos U 



Vl + (r-r) 2 



we can easily compute that dO|cj- ^ and that the line element ds 2 of the metric 

9 = — T 9FRW 

cr 

restricted to 9~ is of Bondi form, namely, 

ds 2 ^- =2d£ldv + dS 2 . 

Here, v = r + r and d§> 2 denotes the standard measure on the two-sphere § 2 , dO 2 + sin 2 6d(p 2 . 
One can check that d T fulfils the conformal Killing equation 

£>d T 9FRW = -2(d T \n(a))g FRW 

and is, hence, a conformal Killing vector of the FRW spacetime. This vector turns out to 
be tangent to 9~ and its restriction to 3" equals ^d v , whose integral curves in turn generate 
Altogether, we find that a conformal Killing vector d T of the bulk FRW spacetime becomes 
a proper Killing vector on its boundary 9~ endowed with the Bondi metric. However, the full 
symmetry group Gq- of 9~ is much larger than that. In the asymptotic power-law case, it 
is given by symmetry group of the conformal boundary of Minkowski spacetime, the so-called 
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Bondi-Metzner-Sachs group (see for instance [DMP05]). In the asymptotic de Sitter case, 
however, the situation is different since the FRW spacetime can be extended to and beyond 
Q~ without a conformal rescaling 1 . In fact, it turns that the symmetry group in this case is 
constituted by the set of diffeomorphisms of R x S 2 which map the point (v, 9, <p) to (e^'^v + 
P(0,(p),R(9,<p)) with a,P € C°°(S 2 ) and R G 50(3) [DMP07]. A direct inspection of this 
formula reveals that the so-called horizon symmetry group is given by the iterated semidirect 
product 

GQ- = 50(3) x (0°°(S 2 ) k 0°°(S 2 )). (5) 
2.2 Classical Scalar Fields on FRW Spacetimes 

This subsection shall offer a brief introduction of both the symplectic space and the mode 
decomposition of a classical scalar field on a FRW spacetime with flat spatial sections. Although 
most, if not all, the material has already appeared elsewhere, we feel that it is worth devoting 
a few lines to this topic as a guideline for the discussion of the more complicated case of Dirac 
fields. 

Hence, let us consider the field : M — > R, whose dynamics are ruled by the conformally 
coupled Klein-Gordon equation 

P<f> = 0, P = -a + ^R + m 2 , (6) 

where R is the scalar curvature, while m 2 > is the squared mass. Since all FRW spacetimes 
as in (1) are globally hyperbolic, one can associate to (6) a meaningful Cauchy problem. Par- 
ticularly, each smooth and compactly supported initial datum yields a solution <p G 0°°(M). 
The linear set of such solutions can be characterised as S(M) = {(f>f \ P(f)f = and 3/ G 
Cg°(M), (f) f = Ef}. Here, E = E~ — E + : C^{M) -»• C°°(M) stands for the causal propaga- 
tor, the difference between the advanced and the retarded fundamental solution. Furthermore, 
for each 4>f G S(M), supp 4>f C J+(supp (/)) U J~(supp (/)), and S(M) forms a symplectic 
space if endowed with the following weakly non-degenerate symplectic form, 

a M {^A 9 ) = j (tf^-^f-h) a(r) 2 d 3 x = E(f,g). 

Notice that om is independent of r since it is in general independent of the Cauchy surface 
constituting the integration domain. 

Although the above characterisation of the space of classical solutions of (6) is, apart from 
the special but inessential choice of Cauchy surface in the definition of ctm> valid on every 
globally hyperbolic spacetime, we can provide a more detailed characterisation of S(M) in the 

1 Note that the conformal transformation relating qfrw and g which made g regular and of Bondi form on 
Ss - is proportional to the identity on in the asymptotic de Sitter case. Nevertheless, we have not omitted it 
in favour of a uniform notation. 
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present case. To wit, by making use of the big isometry group of (M, gFRw)i we can specify each 
element of S(M) via an explicit decomposition in Fourier modes. Particularly, the following two 
propositions have been proven in [DMP07, DMP09]: 

Proposition 2.1. For every k G [0, oo), let T k (r) be a complex solution of the differential 
equation 



{8 Z T + k + a z m z )T k = , (7) 
which fulfils the following normalisation condition 

T k d T T k - T k d T T k = i . 

Suppose, furthermore, that the functions k i— >■ T k (r) and k i— >■ d T T k (r) are both polynomially 
bounded for large k uniformly in r and in L 2 ([0, fc], fcd/c) /or every k > 0. Under these hypotheses, 
every element <fi o/S(M) can 6e decomposed in modes as 



<f>(r, x) = J ^(r, x)<K£0 + ^(r, f)^(fc)J (8) 
where <p^{r,x) is solution of the conformally coupled Klein-Gordon equation of the form 

As shown in [DMP09, PinlO] modes fulfilling the assumptions of the above proposition can 
be concretely constructed in the class of spacetimes we consider via a converging perturbation 
series; moreover, they can be chosen in such a way that the following initial conditions are 
satisfied 

lim e ^T fc (r) = -L, lim e ikT d T T k (r) = -iJ\ . 



The second proposition deals instead with the invertibility of the decomposition (8). 

Proposition 2.2. Under the hypotheses of the preceding proposition, for every 4> £ §(M) ; (8) 
can be inverted as 

4>(k) = ia M {4> k , 4>) ■ 

Furthermore, 4>{k) is square integrable and, as \k\ diverges, it decays faster then any inverse 
power of \k\. 

Notice that, since the symplectic form does not depend on r and since (f)r is a solution of (6), 
4>(k) is time-independent as well. This observation and the preceding discussion entail 



<T M (<t>i,<h) = ~i / Mk)Mk) - Mk)Mk) d 3 k 
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which is once more manifestly time-independent. 

To conclude this section and as a later useful tool, we recall that also the causal propagator 
assumes a somehow more manageable form when decomposed in modes, namely [LuRo90], 

E(T x ,x,T y ,y) = -i J (/>k(T x ,x)(f>k(ry,y) - (j)^T x ,x)(j)^Ty,y) d 3 k . 
2.3 Classical Dirac Fields on FRW Spacetimes 

We shall now focus on the analysis of the classical dynamics of Dirac spinors living on the class of 
FRW spacetimes we consider. The full discussion of the well-posedness of this dynamical system 
requires the introduction of several additional structures and, hence, a rather lengthy detour. To 
avoid such a detour, we point the interested reader either to [DHP09], where an extensive review 
is available, or to the appendix A which contains a small resume of the employed notations and 
conventions. 

As described in the above-mentioned appendix, a Dirac spinor ip (respectively cospinor ip') in 
a given spacetime M is a smooth global section of the Dirac bundle DM (respectively dual Dirac 
bundle DM'). As discussed in lemma 2.1 of [DHP09], DM is well-defined and trivial whenever 
M is globally hyperbolic and simply connected, thus, on a FRW spacetime in particular. In 
other words, in {M,gpRw), the standard picture of (co)spinors being vector-valued functions 
ip : M — > C 4 and ip' : M — > C 4 is valid. That said, we call dynamically allowed any ip which 
fulfils the so-called Dirac equation 

Dip = (-y + ml 4 )ip = 0, 

where I4 stands for the four dimensional identity matrix, whereas y = 7 a e^V^. Here 7 , 
a € {0, 1,2,3}, are 4x4 complex matrices which form a representation of the Clifford algebra 
(7/(1,3) and which are usually referred to as 7— matrices. Due to the employed signature, the 
explicit form of the 7-matrices slightly differs from the standard one; if we denote by o~i the 
standard Pauli matrices, we can set 

^=<o -l)> ^(4 °i) ■ < io > 

Henceforth, lower-case Roman letters mean that the associated quantities are expressed in terms 
of a Lorentz frame e a , a £ {0, 1, 2, 3}. This is a set of four global sections of the tangent bundle 
which fulfil the relation g(e a , e&) = r] a b, where r\ stands for the flat Minkowski metric. All lower- 
case Roman indices are thus raised and lowered via 77^. Furthermore, V denotes the covariant 
derivative on (tensor products of) the full Dirac bundle (see definition 2.10 and lemma 2.2 in 
[DHP09]), its explicit form thus encompasses the spin connection coefficients. 

Remark. As a notational convention, whenever a 7-matrix appears with a definite index, this 
refers to its expression in the non-holonomic basis, e.g. 7 = 7 a | a= o- If we choose the conformal 
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coordinates displayed in (2), then the 7-matrices 7^ are related to the former via a multiplicative 
factor cl(t)~ 1 . 

Let us now recall that, according to theorem 2.3 in [DHP09], the Dirac operator D admits 
unique advanced and retarded fundamental solutions which are continuous linear maps from 
D(DM), the set of compactly supported smooth sections of the Dirac bundle DM into £,(DM), 
the space of smooth sections. Furthermore, it holds for all / G D(DM) that 

DS ± = id D(DM) = S ± D , supp (S ± f) C J^supp (/)). 

Hence, we can introduce the causal propagator S = S~ — S + associated to D, which in turn 
allows us to characterise the space of smooth solutions of the Dirac equation with compactly 
supported initial data as 

6(M) = {tp f I Di(; f = and 3/ G V(DM) such that V/ = Sf} . (11) 

The set &(M) carries a natural Hermitian structure, namely, for all /i,/2 G T>(DM) we can 
define an inner product as 

= i J 4vt*hdiA(V) = -iJ flSfr dfi(M), 

S M 

where £ is any Cauchy surface of the underlying background whose unit normal vector is denoted 
by n. Such choice plays no role since independence of the integral from S was already proven 
in proposition 2.2 of [Dim82]. Here, the superscript f refers to the Dirac conjugation map 
which, for all ip G £(DM) is defined as ip^ = if)* j3 where * denotes adjoint with respect to the 
standard inner product on C 4 , whereas f3 is the Dirac conjugation matrix, that is, the unique 
element of SX(4, C) which fulfils 

(3* = 13, 7* = -/37a/?- 1 , 

and — i/3n M 7 M is a positive definite matrix for all timelike and future-pointing vector fields n. 
In the representation (10) we have chosen, j3 = —ijo- 

Remark. We only concentrate on Dirac spinors since the behaviour of cospinors can be inferred 
from the one of spinors by applying the Dirac conjugation map to all relevant equations. 

We now specialise our treatment to any FRW spacetime {M,g FRW \ with flat spatial sections 
whose scale factor a(r) is of the form (4) . To analyse the behaviour of a Dirac field living thereon, 
we first discuss the choice of a Lorentz frame, and two natural possibilities exist. On the one 
hand, since the spacetime (M, gFRw) can be conformally embedded 2 in a larger spacetime 
(Me, gE) which contains the horizon we could simply choose the frame e a as the global one 

2 In the asymptotically de Sitter case, the FRW spacetime at hand can even be isometrically embedded in a 
larger space which, however, differs from the Einstein static universe (Me^e)- 
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on Me- The advantage of this choice would be that this Lorentz frame would be automatically 
well defined on On the other hand, in the following discussion it will be crucial to exploit 
the conformal flatness of qfrw displayed in (2). Hence, we shall choose the standard Minkowski 
frame, say e a , and then define the one on (M, gFR\v) as 

e a = a(r) _1 e a . 

Notice that, due to the divergent factor a, the frame e a cannot neither be extended to 9~ nor 
to (Me, 9e)- Although we will be able to circumvent this problem implicitly, we refer the reader 
interested in further explicit details to [HalO]. In the chosen conformally flat Lorentz frame, 
the covariant derivative on DM, and, consequently, the Dirac equation assume a somehow more 
manageable form, viz. 

a~l (-7°<9 T - YVi + am) a? ip = , (12) 

where Vj is the Cartesian gradient along the spatial directions. 

As in the scalar case, we would like to have a more explicit Fourier mode characterisation 
of solutions of the Dirac equation on (M, gFRw) with compactly supported smooth initial data. 
To this avail, and in view of the high symmetries of a FRW spacetime, it would be desirable 
to reduce (12) to a diagonal form. The natural approach which calls for diagonalising D'D is 
somehow not well-suited to our aims, since, after performing the Fourier transform along the 
spatial directions, the resulting ODE displays a term which is both complex and linear in k, a 
situation which leads to potential practical difficulties when constructing a state. A different 
and ultimately more effective diagonalisation procedure has been introduced in [BaDu87] and 
we shall now discuss an improved and more clear version of this approach. To start, let us 
introduce the modified operator 

^ ■ On -M d T + iam -V ■ a \ a 
\ —V ■ a o T — tarn J 

where Oi are the Pauli matrices and • denotes the inner product in M 3 . The equations Dip = 
and Dip = are equivalent, as — 7 is invertible. If we define D' = — aD'7 , then we can compute 

^ " a \ (-0? + A + ia'm + a 2 m 2 )I 2 ) ' (i4j 

where / denotes a derivative with respect to r and A is the Laplace operator on 1R 3 . Notice 
that, on Minkowski spacetime, the above procedure is equivalent to considering DD' , which is 
already diagonal in this special case. The following lemma summarises the above discussion and 
will play a key role in the study of the fall-off behaviour at infinity of solutions of the Dirac 
equation: 
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Lemma 2.1. The following Cauchy problems for a spinor tp : I x R 3 — > C 4 are equivalent: 

' DD'u = (- V M W + | + m 2 )u = 
2) ^ u\ t=0 = 

k (9f«)|t=o = 7o/ 



1) 



= o 



3) 



3)^ = 







4) 



Pjjtt = (J)S)')u = 
n| t=0 = 



t=o 



/ 

a[0 



where ip = D'u and ip = D'u. 

Proof. 1) <^=^ 2) has been already proven in theorem 2.2 of [DHP09]. 1) <^=^ 3) is immediate 
once we recall that D differs from D by a multiplicative pre-factor which is an element of 
SL(2, C°°(R, (0, oo))). 3) <^=^ 4) can be proven as follows: suppose 4) holds true, and let us 
introduce tp = ID'u, then ^(0) = D'u(0) = 0(0)70-0512- The restriction of u and of its derivatives 
to S is meaningful since is, per direct inspection of (14), a set of four second order hyperbolic 
differential operators. Thus, u is a unique smooth solution on the whole M. In order to achieve 
compatibility with the initial condition of 3), one sets ^(0) = / and thus a(0)7o7 M V At : u|s = f ■ 
Since V^njs = —n^^ where n is the normal vector field to the Cauchy surface, then, owing to 
n = dt, one obtains per direct substitution — a(0)7o7°^-|£ = /, which, thanks to the assigned 
initial condition for 4) and to the identity 7o7° = 1 yields the sought result. Notice that we 
have just proved that 4) implies 3), but uniqueness of the solution of the Cauchy problem also 
entails the converse. The transitivity property of the proven equivalences suffices to infer the 
statement of the lemma. □ 

The main advantage of this lemma is that, in order to construct a solution of the Dirac equation 
in a cosmological spacetime, we can just focus on the fourth Cauchy problem. Furthermore, 
using the fact that the operator DD' is diagonal, we can discuss the four components of the 
solutions u independently. It turns out that a convenient mode basis of solutions of DD' is given 
by 



UkA T ) e 



(2vra): 



(15) 



where k = \k\ and 

/ Xk,i 


u M = Q 

V o 

/ Xk,2 \ 



V o / 



/ \ 



Uk,2 





Xk,i 





Uk,3 



Uk,6 



V o / 
/ \ 

Xk,2 



V o / 



Uk,7 



( o \ 



V o / 



/ \ 





UkA 



Xk,2 

\ o / 



Uk,S 



( o \ 





V XM / 

o \ 



\Xk^ J 



( 
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Here, Xk,i an d Xk,2 constitute two linearly independent solutions of the ODE 
?Xkj = (dr + k 2 + a 2 m 2 - ia'm) X k,j(j) = . j £ {1,2} 



(16) 



Remark. Although 7 is not a real differential operator and, hence, Xk,i an d Xk,2 cannot be 
related by complex conjugation, we are still free to choose XkJ an d Xk^2 as a basis of the space 
of solutions of 3 . 

As in the study of the solutions of the Klein-Gordon equation conformally coupled to scalar 
curvature, we would like to use (15) to construct a mode decomposition of each solution ip of 
the Dirac equation. Yet, in view of lemma 2.1, which guarantees that a solution of the Dirac 
equation can be found from one of DD' by application of 2)', we can expect that not all 8 basis 
modes provided by (15) are needed to construct a mode expansion of ip. As a matter of fact, the 
following result implicitly entails that already the first four modes listed in (15) are complete. 

Proposition 2.3. Let X k be a smooth solution of (16) which satisfies the following normalisation 
condition 

\(d T + iam) Xk \ 2 + k 2 \ Xk \ 2 = l. (17) 

Suppose, furthermore, that for any fixed r and k\ > 0, both the functions k >->■ X k(r) and 
k i y d T Xk{T) are in L 2 ([0, k±], k 2 dk) and they grow at most polynomially in k. Then, for any 
solution ip £ &{M), it holds 

4 r~ 

V(r,f) = Y J J Mk)^(T,x)d 3 k, (18) 

where 

^ kl (T,x) = V> U ) l{T)e 7 X . (19) 
(2vra)2 

If one takes into account (17), it is also possible to invert the above decomposition as 

Proposition 2.4. Under the hypotheses of the preceding proposition, the mode decomposition 
of any Dirac field if) £ &{M) can be inverted as 

i>i{k) = -i j a 3 (r) V^(r, x) 7o V(r, x) d 3 x = s M (^,v ^) ■ ( 20 ) 

Furthermore, each ipi(k) is a square integrable function which decays faster then any inverse 
power of k. 
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Proof. In order to show that either (20) inserted in (18) or vice versa yields an identity is a 
matter of a long, tedious, but direct computation. However, the statement on the regularity 
of ipi(k) requires a closer look. It descends both from the hypotheses on Xk formulated in the 
previous proposition and from the fact that, at fixed r, ip(r,x) is an x-dependant four-vector 
whose components are compactly supported smooth functions. Hence, they are square integrable 
and their Fourier transform is rapidly decreasing in k. If one combines this result with the fact 
that ipi(k) is a linear combination of the product of the Fourier transform of compactly supported 
smooth functions together with either imaxk(T~), (cr ■ k)xk(T~) or d T Xk(i~), then the sought result 
follows. □ 

Notice that, as in the case of a scalar field, the r— independence of ipi(k) stems from both the 
Cauchy surface independence of 5m and the fact that ipr , is a solution of the Dirac equation. 

Remark. For the class of spacetimes under investigation, a set of solutions of (16) Xk, which 
fulfil the properties required in the above two propositions, are constructed in lemma B.l. They 
are uniquely determined by the initial conditions we shall assume to be fulfilled in the following, 

viz. 

lim e ^ Xfc ( r ) = _L l im e ikT d TX k(r) = ~4= • 



T— > — OO 



V2k' r--^ ^2 



As a first profitable consequence of the mode decomposition we can obtain a convenient expres- 
sion for the conserved Dirac Hermitian product on FRW spacetimes with flat spatial sections. 

Lemma 2.2. For every ipi and ip2 in &{M), it holds 

4 ~ 

e M (il>i,ik) = 52 J Mk) ^, )2 (fc) d 3 k, (21) 

1=1 R3 

where $ij(k) =s M {^ v ^j) for j G {1,2}. 
Proof. Let us rewrite 

SmWi,^) = -i j iplrffadnCE) = J iplfa a 3 (r) d 3 x 

where we have used both ijv = ip* (3 and the freedom to choose S in such a way that n = dt and 
— i(3rf = I4. We now insert the expansion 



ip 2 (T x ,x)= / dk^s M {i>%^^2) ip^(T x ,x) 

to>:>, '=1 



into Sm^i,^)- Since the x-integration is over a compact set while the integrand is rapidly 
decreasing for large k, we can use the theorem of dominated convergence to switch the order of 
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the k- and x-integration. Afterwards, by using the sesquilinearity of Sm and by noticing that 



To conclude the section, we remark that the following mode decomposition of the causal 
propagator can be obtained: Provided that the modes Xk fulfil the conditions imposed in the 
two main propositions of this section, one finds [HalO] 



In the last section of this paper we will use this expression to derive a mode expansion for the 
two-point function of the states we shall introduce. 

3 Quantum Field Theory in the Bulk and on the Cosmological 
Boundary 

The aim of this section is to use the constructions of the previous section in order to achieve two 
main results. On the one hand, we will show that it is possible to set up a genuine Fermionic 
quantum field theory on past null infinity 9f - . On the other hand, we shall prove that one can 
encode the information of a bulk quantum field theory into the counterpart on 

As in the previous section, we shall first discuss the case of a massive real scalar field confor- 
mally coupled to scalar curvature in order to recapitulate the main features of the construction 
we aim for in a scenario where it is already known to work. Subsequently, we prove that the 
same procedure can be successfully carried out in the case of a Dirac field as well. 

3.1 The Bulk-to-Boundary Correspondence for Real Scalar Fields 

This subsection recollects some of the results already proven in [DMP07, DMP09, PinlO] con- 
cerning the projection of real scalar fields to The key step is the following: while the 
classical theory on (M, gFRw) naturally contains a symplectic space of solutions of the equa- 
tion of motion, there is no dynamical content on the cosmological horizon. Therefore, we are 
forced to perform a choice of a symplectic space on the boundary which can be only justified a 
posteriori. To wit, as in the above-mentioned references, we introduce 



where d§> 2 (6,ip) is the standard measure on the 2-sphere. This is a strongly non degenerate 
symplectic space if endowed with the symplectic form 



Sm{4>i, V'jfc/) = , the sought result follows. 



□ 




S(9T) = {$ G C°°(3T) I $ and G L 2 (^-;dvdS 2 (9, tp))} , 




V$i,$ 2 G S(9T) 



MxS 2 



16 



The motivation to choose (S(9f ),<Tc s -) is the existence of a symplectomorphism 

T s : S(M) -»• S(Qf-) 
whose explicit action on the elements of S(M) is defined as 

Fs((p)(v , 9, if) = lim —ua(j)(T(u,v),x(u,v,0,(p)). 

u— >— oo 

Here, we have first implicitly switched from the Cartesian coordinates x = (x, y, z) on the bulk 
Cauchy surface isomorphic to M 3 to the spherical ones (r = \x\,9,ip) and introduced the null 
coordinates v = r + r and u = t — r afterwards. The above limit was already computed in the 
proof of theorem 4.4 in [DMP07] for the asymptotic de Sitter case, but the same computation 
can be trivially repeated for the asymptotic power-law case, as it relies on the fall-off behaviour 
of the classical modes towards 9~ which is the better the larger 5 in (4) is. The result found in 
[DMP07] is 

T s (<f>)(v, 0,<p) = -j= J [e~ ikv cf>(k, vr - 9, vr + + e +ikv cf>(k, vr - 9, ir + ^)) - dk . 

Remark. The ansatz for the map T,s is based on the following preliminary considerations. 
In Minkowski spacetime, smooth solutions 4> of the massless Klein-Gordon equation are known 
to decay as (— u)^ 1 towards that is, in the limit u — > — oo. To wit, due to the fact that 
Minkowski spacetime can be conformally embedded into the Einstein static universe and that 
the scalar field has conformal weight 1, 

2 

6 = : <b 

^J(l + u 2 )(l + v 2 ) 

is a solution of the massless, conformally coupled Klein-Gordon equation in Me [DMP05]. 
Hence, <f> is manifestly finite on Q~ and the decay behaviour of <p follows. The effect of the 
projection map T$ can now be understood as follows. The multiplication with the conformal 
factor a transforms the massive, conformally coupled field <j) on (M, gFRw) m t° a scalar field on 
Minkowski spacetime with mass am. This mass term is, however, not essential for the fall-off 
behaviour of the associated solutions, as it is finite or vanishing on G - . Hence, the additional 
factor of — u in Ts is both necessary and sufficient to obtain a finite projection of solutions in 
(M,g FRW ) to 3r. 

Since 4> decays faster than any inverse power of k and since it is square integrable, Tg(<p) is 
an element of S(9~). Furthermore, as anticipated above, the following proposition (see theorem 
4.2 in [DMP07]) holds 

Proposition 3.1. The symplectic form om is conserved under the projection Ts, that is, for 
every 4>\ and 4>2 in S(M) ; it holds 

ffM^i,^) = o-^-{T s 4>i,T S (t)2) ■ 
Thus, the map T is an injective symplectomorphism from S(M) into S(9~). 
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This proposition relating classical field theories entails a close relationship between the bulk and 
the boundary algebra of quantum observables. As a matter of fact, let us recall that to every 
real vector space S endowed with a weakly non-degenerate symplectic form as one can associate 
a unique (up to isometric *-isomorphisms) Weyl C*-algebra whose generators W(s), s G S, fulfil 
the defining relations: 

a)W(s)* = W(-s) b)W{s)W{s') = e^ asi - s ' sl) W{s + s'), Vs,s'g5. 

In the scenario at hand, it is thus clear that it is possible to associate a Weyl C*-algebra both 
to the boundary - W(9~) - and to the bulk, say W(M). Furthermore, in view of the above 
proposition it holds (see theorem 4.2 in [DMP07]) that 

Proposition 3.2. The symplectomorphism T$ induces an injective *—homomorphism of C* — algebras 
ts : W(M) — > W(Q~), fully determined by 

i s W{<p) = W{T S (p) , G S(M) . 

Furthermore, for every state, that is, for every continuous, positive, normalised, linear functional 
co : W(G _ ) — > C, there is a counterpart ojm '■ W(M) — > C unambiguously defined as 

co M {W {<!>)) = co {is{W{4>)) = u {W(T s <f>)) • 

3.2 The Bulk-to-Boundary Correspondence for Dirac Fields 

In this subsection we will show how the bulk-to-boundary correspondence introduced for the 
scalar field can be adapted to Dirac fields. To this avail, the first tool we need is the counterpart 
of the Dirac bundle on the boundary. Although Q~ could be considered as a genuine manifold 
on its own, the question whether it admits a spin structure can not be answered automatically 
with the results we have invoked up to now, as 9~ endowed with the Bondi metric is not 
a globally hyperbolic spacetime. However, if we recall that 9~ is defined as a codimension 
1 embedded submanifold of the globally hyperbolic Einstein static universe (Afg,^), we can 
exploit definition 5.3 in [Hus96] and introduce = l*(DMe), where i : 9~ Me is 

the natural embedding map. Furthermore, on account of corollary 6.7 in [Hus96] and of the 
isomorphism 3 DMe ~ Me x C 4 , we can conclude that Z?Q~ ~ Q~ x C 4 ; consequently, the 
associated smooth sections can be equivalently interpreted as elements of C°°(3 ; ~,C 4 ). 

As in the previous subsection, we have to specify a suitable boundary configuration space 
whose choice is justified only a posteriori. To wit, we call 

©(9T) = G C°°(Qr, C 4 ) n L 2 (D%-;dv dS 2 ) , H(a)*(u, n) = V(v, n)}, (22) 

where L 2 (D9~; dv <i§ 2 ) is isometric to L 2 (Q - ; dv dS 2 ) <8> C 4 . Furthermore n G M 3 has unit 
Euclidean- norm and thus corresponds to a point on the two-sphere, whereas H(n) is the following 
matrix-operator 

^ (A R *) ■ < 23 > 

3 Note that DMe is trivial since (Me,9e) is both globally hyperbolic and simply connected. 
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where n-a is the Euclidean inner product of n and the Pauli vector. Notice that, in sharp contrast 
to the scalar case, a reminiscence of the bulk equation of motion prevails on the boundary via the 
constraint 3(n)^f(v,n) = ^/(v,n). This identity implies that only two of the four components of 
an element of ©(3~) are independent. 

©(3? _ ) can be naturally endowed with the Hermitian form 



Rx§ 2 



where, again, vj/t = — i\l/*7o while * is the standard adjoint map on C 4 . This Hermitian form is 
non-degenerate on ©(Q 5- ) because (n • a) 2 = I2 and, hence, 

a 9 -(*i,* 2 ) = 2 J (^1^2)1 dvdS 2 , 

RxS 2 ie { 1 ' 2 > 

where only the two independent components of \£* and ^2 have been summed. Furthermore, 
since 9~ is topologically equivalent to R x S 2 , we are free to Fourier-transform each G ©(S~) 
as 



I €{1,2,3,4} 



Accordingly, the Hermitian product reads 



r 4 ^ _ 

s 9 -(*i,* 2 )= J ^*i,i(k,0,<P)*iMW,v)dkdS? '. 

RxS 2 1=1 

The main result of this subsection is the following proposition which justifies the choice for 
©(S - ). Regarding coordinates we employ the same notation as in the previous subsection. 

Proposition 3.3. The map T D : S(M) -»• ©(Of - ) 

FdC^) = lim —ua 3 / 2, ip(T(u,v),x(u,v,6,tp)) 

U — > — OO 

possess the following properties: 

a) Td is a well-defined, that is, Fd(<3(M)) C ©(3~), 

&,) //we se£ ^ = r^)^) and denote the four-vector representing its Fourier transform by ^, 
it holds 



tf(M) = \k\ (I4 + H(n)) 





(ijji (k,-n)\ 




( ° ^ 




&(k) 


1P2 (k, -n) 



+ ®{-k) 




Tp 3 (-k,n) 




\ 


^ ) 




\ij>4 {-k,n)j 


) 



(24) 



where ft is still a unit vector in R 3 , S is the matrix defined in (23), and ipi(k) are the mode 
expansion coefficients of (18), 
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c) preserves the Hermitian forms, that is, 

SM(ipi,ih) = s 3 -(r D (^i),r D (^ 2 )) 



Proof. To prove a), we follow once more the proof of theorem 4.4 in [DMP07]. We start from 
the identity 



hm (-u)a 3 ^(u,v,e,ip) = lim — — / V i)'u ktl (r)e ik %(k) d 3 k , 

u — y — oo u—i-—oo (ZTT)''' J z — ' 

3 3 — * a A 

where ID' = a^T)'a~^ . If we write k in spherical coordinates k,9,ip where 9 measures the angle 
between x and k and if we restore the radial coordinate so that v — u = 2r, we obtain 

9(v,0,<p) = lim [ i)' Uk:l (r)e ikrcos ^Mk)k 2 sm(0)dkd9d^. 

R+x§ 2 

On account of the properties and estimates of \k discussed in 2.3 and appendix B and, 
setting c = cos(0), we obtain 

2-7T 1 OO 

*(v,9,v) = Jjm^—^J J J (p' s {k)e s ^e ikrc + 0(\uD) $ s (k) ^ dk dcdtp, 

o-io *={+-} 

where ip + = (V>i, 0, 0), ip- = (0,0,^3,-04) and ip = ip + + tp-. In the preceding expression 
D'±(k) is the 4x4 matrix 

, (f , (25 ) 

while the 0{\u\~ e ) contribution descends from the estimates in appendix B and, in the asymptotic 
de Sitter case, from an expansion of xt( T ) v &hd for large r that can be derived from equation 
8.421(9) in [GR95]. In the asymptotic power-law case, the 0(|u|~ e ) contribution is identically 
vanishing. 

The multiplicative u factor in front of the limits can be cancelled via an integration by parts 
in c. Afterwards, also the term proportional to 0(\u\~ € ) vanishes in the limit of divergent u 
thanks to dominated convergence. What is left after these operations consists of three terms, 
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each to be evaluated as u — > — oo: 

27T OO 







2n oo 

— It 



l ~v) (2^// S ^(fc)e a( - ifcT) e ifcr ^(A;,0^)dA;#+ 



2r(u. „ . 

«€{+.-} 



2ir 1 so 



J J J Y D' s (ky { ~ ikT) e ikrc 9Mk,c,<p) dkdcdtp 



2r(u, v) (2tt) 3 / 2 „ 

o-io «e{+,-} 

The second and the third term can be seen to vanish in the large u-limit by an application of the 
Riemann-Lebesgue theorem since the ratio —u/(2r) tends to 1 and d c tpi is a regular integrable 
function. Notice that, in the last term, there is a potential obstruction in the application of the 
Riemann-Lebesgue theorem due to the behaviour of the integrand at c = —1. Yet, this problem 
can be overcome if one splits the domain of integration in [— 1, — 1 + e] and [— 1 + e, 1]. The 
second contribution vanishes for large — u and for all e > 0, while the first one yields a finite 
result whose value is regulated by e itself and, hence, vanishing as e — > 0. We are left with the 
first of the three integrals and, since both ip s (k,O,0) and ip s (k,n,(p) do not depend on <p, this 
variable can be integrated out, yielding 

oo 

*M,¥>) = [ dk D' + (-k)e- ikv ^+(-k) - D'_{k)e ikv $-(k), 

V27T J 


where k = (k, 6, ip). Due to the form of D±, the previous expression can be rewritten as 

oo 

n) = (I 4 + S(n))-^= j dk k [e~ ikv ^+ (k, -n) + e ikv ^ (k, n)) , 
o 

where S(n) is defined in (23). Thanks to the regularity of ip stated in proposition 2.4 and 
regardless of the multiplication by k, every component of ^ is a smooth square-integrable func- 
tion. This allows to conclude the proof of a) by noticing that the constraint in the definition 
of <5(Q~) is fulfilled by since H(n)(l4 + H(n)) =14 + H(n). Statement b) holds because the 
Fourier transform of the previous expression becomes equal to (24) if we change the variable k 
in —A; in the second contribution to the integrand and introduce the appropriate Heaviside step 
functions. 

In order to prove c), we insert the Fourier decomposition of r^^i) and ^(^2) found in b) 
in the definition of s$- to obtain 

r 4 ~ 

s 9 - (r D (v>i), r D {*h)) = I Y ^M fc . e > <p)MK 0i v)k 2 dk d ® 2 ■ 

R+x§ 2 1=1 
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The right hand side of the above equality coincides with %(^i,^2) evaluated in Fourier space 
and in spherical coordinates. □ 



Remark. In this case, the ansatz for is motivated as follows. As the conformal weight of a 
Dirac field is |, the factor a 2 in transforms the massive Dirac field in (M, gFRw) a Dirac 
field in Minkowski spacetime with mass am. The fall-off behaviour of corresponding solutions 
with compactly supported initial data is again independent of this time-dependent mass term. 
Since the massless Dirac field in Minkowski spacetime can be understood as a collection of four 
massless scalar fields, the factor — u in To is, as already discussed in the context of F$, necessary 
and sufficient to obtain a well-defined projection of solutions to 

The last proposition is the analogue of proposition 3.1 for the Dirac fields, namely, the key 
building block to establish a relation between the bulk and the boundary algebra of quantum 
observables. Hence, let us start introducing these algebras following the discussion presented in 
[Ar70] and further developed in [Sa08, DHP09, SanlO, HalO]. The essential idea is that, in order 
to describe the field algebra of Dirac fields 'B(M), it is convenient to treat spinor and cospinor 
fields as a single and combined object. We refer to [Sa08, DHP09] for an in-depth analysis of 
the construction and analysis of this double-algebra and only summarise the key steps in the 
following. To wit, by ¥>(M) we indicate the field *— algebra generated by the identity and the 
linear functional B(f) with / = /1 © f 2 G T>(DM © DM*) subject to 

B(Df®D'h) = 0, B*(f) = B(Tf), 

where Tf = (f 2 ,fl) with f denoting the Dirac conjugation. Furthermore, all generators are 
required to satisfy the anticommutation relations 

{B(f),B(h)} = iS®(f,h) 

where S®(fi © f 2 , hi © h 2 ) = iS(h 2 , fi) + iS(f 2 , hi) = s M (V'/i . ^h 2 ) + *M (ipf 2 , ) and where S 
is the causal propagator of the Dirac equation. In this framework, we can recover spinors and 
cospinors as ^(/l) = B (f) with / = fi © and ip(f 2 ) = B(f) with / = © f 2 . 

The assignment of a field algebra to the theory on the horizon is slightly more involved. Let 
us start noticing that, if we define ^ on as — i^*7o for all ^ G 6(9"), we can mimic the 
construction of 'B(M) to construct its boundary counterpart H>(^s~). To wit, following [HalO], 
we define 6®(9~) = 6(9") © 6(9~)t where 6(G~) t is the set formed by ^ for every ^ in 
(5(9~). can be equipped with the inner product 

s|_ ((ai,4),(>i,4)) =s Q -{b 2 ,ai) +Sc ! -(a 2 ,b 1 ) 

and the conjugation j 

j(ai,4) = (a 2 ,a|) . 
We summarise the above discussion in the following definition. 
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Definition 3.1. The boundary algebra of Dime fields 23(9 ) is the topological *— algebra gener- 
ated by the linear Junctionals f i-> B(f) with f G (5®(9 _ ) together with the following conditions. 

a) The fields fulfil the canonical anticommutation relations (CAR) 

{B(f), B(g)} = 4_ (/, g) V/, g € 6® (Or) . 

6j TTie *-operation is specified by the antilinear involution * : 23(Q~) — > 25(9") defined as 

B(f)* = B(j(f)) V/€6®(9T). 

23 (9~) is endowed with the quotient topology which descends from the local Frechet topology 

°/e n (c°°(3- c 4 ))®». 

On the horizon, there is no equation of motion to be taken into account, but a reminiscence of 
the Dirac equation has been implemented in the construction of 6®(9~). 

We are now ready to formulate the main proposition of this subsection which establishes a 
relation between the boundary and the bulk algebra. The proof is a straightforward application 
of the previous definitions and of proposition 3.3. 

Proposition 3.4. The map id '■ 23(M) — > 23(Q~) unambiguously determined by its action on 
the fields B(f), f G T)(DM © D*M) as 

i D {B(h © /])) = B(r D y> fl ) © r D % 2 )t) V / = h © / 2 e t>(dm © d*m) 

is an infective *—homomorphism. Furthermore, for every state, that is, for every continuous, 
positive, normalised, linear functional u : 23(9 _ ) — > C, there is a counterpart ojm '■ 23(M) — > C 
unambiguously defined as 

uo M (B{f)) = wM£(/i © / 2 )) V/ = h © h G ?>(DM) © V(D*M) . 

4 Boundary States and their Hadamard Bulk Counterpart 

In the previous section, we have proven that it is possible to induce a state for a scalar and 
Dirac quantum field theory on the bulk spacetime (M, gpRw) by a state defined on the relevant 
quantum theory living on past null infinity. However, this construction of a distinguished state 
would be moot if the result could not be shown to be physically meaningful. Particularly, it is 
mandatory to check if the constructed state, say com, fulfils the so-called Hadamard condition. 
From a physical point of view, this condition guarantees that the UV behaviour of lom mimics 
the one of the Poincare-invariant Minkowski vacuum and, hence, that the quantum fluctuations 
of observables, such as the smeared components of the stress-energy tensor, are bounded. From 
a formal point of view, the requirement for ojm to be Hadamard boils down to a constraint 
on the form of the wave front set of the integral kernel proper to the truncated two-point 
function associated of the state ojm [Rad96a, Rad96b, KraOO, Hol99, SaVeOl, Sa08, Sa09]. Since 
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a cohesive and fully comprehensible explanation of all the mathematical tools underlying these 
notions would require a lengthy appendix, we deem it more appropriate to point an interested 
reader to the just mentioned references. 

Our goal is to follow the same procedure employed in [DMP05, DMP07, DMP09, PinlO]: we 
will show, that, mostly thanks to the huge symmetry displayed by past null infinity it is 
possible to assign a distinguished state lo^- to the field algebra on 9" whose pull-back to the 
bulk algebra is of Hadamard form. Furthermore, we will show that, both in the scalar and in 
the Dirac case, wg- can be chosen to satisfy an exact KMS condition on the horizon, a property 
which will turn out to be preserved in the bulk under suitable circumstances. All states we 
construct are automatically quasi-free, that is, their structure is completely determined by their 
two-point functions. However, as recently proven in [Sa09], this requirement is irrelevant for the 
discussion of the Hadamard condition. 

In the case of a scalar field, the assignment of a two-point function ui2 is tantamount to 
providing a distribution on D(M x M). However, in discussing (charge-conjugation invari- 
ant) Hadamard states for Dirac fields, we are concerned with two non-vanishing distributions 
tv + {f,h) = co(tp(h)^ (f)) and tt)-(f,h) = co{^ (f)4>{h)), where / G D(DM) and h G D(D*M). 
In other words, both tu + and ro~ lie in T>'(DM M D*M) and they are related by the anticom- 
mutation relations 

»"(/, h) + to + (f,h) = iS(h,f) 

A quasifree state on the boundary algebra 23 is analogously determined by the choice of 
two distributions which are related by the boundary CAR. 

4.1 The Scalar Field Case 

This case is certainly not a novel one as the quest to construct a bulk Hadamard state by means 
of a boundary one was already pursued in [DMP07, DMP09, PinlO]. Yet, in these papers, 
the possibility of defining a state on 3~ which fulfils a KMS condition was neither mentioned, 
nor analysed. Hence, we shall first briefly recollect the construction in the case of vanishing 
temperature T = (3~ l = case and afterwards extend it to j3 < oo. 

To wit, we adopt here the definition of a quasi-free state introduced in [KaWa91] and consider 
ijj : W(^s~) — > C unambiguously determined by 

w(W($)) = e' ti ^ 1 , V$€S(9T) (26) 

where 

$') = J 2k@(k)$(k, 6, <p)$'(k, e, if) dk dS 2 , 

RxS 2 

and Q(k) denotes the Heaviside step function. Furthermore, $ stands for the Fourier-Plancherel 
transform (see appendix C in [Mo06]) defined for all G S(9~) as 

$(fc,0 )¥ O = --Ly e lkv *(v,9,<p) dv. 

K 
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Notice that \x satisfies the constraint that every quasifree state (see [KaWa91], [M0O6, app. A]) 
has to fulfil, namely 

$01 <4|M*,*)I*IM*', *')!*■ (27) 

One of the most striking properties of u descend from the study of its interplay with the boundary 
group of isometries namely, the BMS group in the asymptotic power-law case, or the 

horizon symmetry group (5) in the asymptotic de Sitter case. In fact, it turns out that u is 
invariant under the *-automorphisms a s : W(Q~) — > W(Q~) on the boundary Weyl algebra 
induced by every element g G Gq- via a fl (W($)) = W^oq -1 ) for all <3? G S(9f - ). Particularly, 
summarising the content of theorem 4.1 in [DMP07], it turns out that u is the unique pure and 
quasifree state on W(Q~) which is invariant under the automorphic action of the one-parameter 
subgroup of Gc-- generated via the exponential map by d v , the generator of the rigid translations 
on Q~ along the M-direction. Moreover, one finds that w is a ground state with respect to the 
^-translations. 

However, as found in [DMP10], the simple construction of the boundary ground state on the 
level of Fourier modes allows to define a boundary KMS state with respect to the translations 
along v by means of the inner product 

00 

/2k — -~ 
l- e -Pk $i(fc,g,y)$2(M,y) dkdS 2 (28) 

—00 

for any /3 > 0. In more detail, the following statements can be proven. 

Proposition 4.1. The two-point function (ip defined in (28) induces a quasifree state up : 
W(Q ) —7- C via up(W(&)) = e 2 which enjoys the following properties: 

a) For every f3 > 0, up is the unique pure quasifree state which is a KMS state with respect 
to v— translations. 

b) up is invariant under the automorphic action of the subgroup Gq- generated by rotations 
and v— translations on W(Q~). 

c) In the limit f3 — > 00, [ip converges weakly to the two-point function of (26). 

Proof. As \fi(&, &)\ $ $)| for all j3 > and for all $ G S(M), up is a well-defined quasifree 

(mixed) state. Invariance under the translations generated by d v can be either directly inferred 
from the explicit form of the two-point function or deduced by repeating the proof of theorem 
4.1 in [DMP07], which additionally entails the uniqueness property. Moreover, invariance under 
the automorphic action of the rotations and v— translations arises from the same reasoning used 
in section 4.2 in [DMP07] . 

Since the state is quasifree, it is sufficient to verify the KMS condition at the level of two-point 
functions. Let us consider G S(9 _ ), and let us introduce the two functions 

F(t) = fipi®, a t (*')) , G(t) = nt,(at($), , 
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where a t (Q)(v, 9, ip) = &(v — t,9,(p). A direct computation shows that both F(t) and G(t) can 
be seen as the Fourier-Plancherel transform of suitable square- integr able functions F(E) and 
G(E) which decay faster then any inverse power of E. Thus they are continuous functions and, 
on account of the explicit form of hb (28), they are related as 

F(E) = e pE G{E) . 

This in turn entails F(t + if3) = G(t) and, hence, the validity of the KMS condition as shown in 
section 5.3 of [BrRo96v2]. □ 

As discussed in proposition 3.2, each boundary state induces a bulk counterpart, particularly, 
we can define 

M ■ 

where is is the injective *-homomorphism introduced in proposition 3.2. As already anticipated, 
we shall now show that u^ 1 fulfils the Hadamard condition as formulated in definition 3.3 of 
[DHP09], which is valid both for scalar [Rad96a] and for Dirac fields [Hol99, KraOO, SaVeOl]. 

Theorem 4.1. The state ujf : W(M) -»• C is a Hadamard state, that is, the integral kernel of 
its two-point function 

is a distribution in D'(M x M) which enjoys the Hadamard property. 

Proof. Per direct inspection, one can realise that the limit (3 — > oo of (28) is the two-point 
function of the state (26) which is of Hadamard form as already proven in [DMP09, PinlO] 4 . 
Thus, in order to conclude the proof, it is sufficient to show that 

h) = uppif, h) - Woo,2(/, h), 

is an element of D'(M x M) with smooth integral kernel. Aiming to prove the continuity of 
first, we set <£>/ = T s {E(f)) and <5> h = T s (E(h)) G S(3~), and directly compute 

oo 
— oo 

where Dp = fjLp — fioo is the difference of the relevant two-point functions of the boundary states. 
Since - is bounded, \Ap(f, h)\ is dominated by the the L 2 -norms of $>f and &h- Using the 
continuity of the causal propagator, seen as a map from D(M) to £(M), the regularity property 
of the modes T^, the continuity of the Fourier transform in L 2 and, finally, the definition of Fs, 
one obtains that 

q(n) 

\\k n $ f \\ L 2<C(K 1 n)Y,\\D 3 fh^, 

3=0 

4 The arguments used in those papers are based on an early proof presented in [Mo08] . 
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where ^ is the Fourier-Plancherel transform on S - , C(K,n) is a constant which depends both 
on the compact set K C M which contains the support of / and on the order n, and H-D^/H^oo 
has to be understood as the L°°-norm of the j-th partial directional derivative evaluated on the 
j directions which give the maximal result in an arbitrary but fixed local coordinate system. 
This observation together with the previous ones yields the sought continuity, viz. 

\^UM <Ck^PV||oo^||^IU, 

jr=0 1=0 

where, once more, the constant Ck depends on the compact set K which contains the support 
of both / and h, while q is a fixed constant. 

Let us now tackle the problem of proving the smoothness of the integral kernel of ojp^ defined 
as the composition Dp o {Tg{E) <g) Ts(E)). To this avail, we can make use of theorem 8.2.13 
in [Ho90] which allows us to control the wave front set of this composite linear functional. To 
wit, let us recall that, in Dp, Jffl 1 decays faster then any inverse power of \k\. Hence, if 
(xi, X2, ki, hi) is a point in the wavefront set of Dp, k\ and k^ must have vanishing components 
along causal directions, i.e. the v-direction, but the restriction of Tg{E) Ts(E) to 9~ x Q~ 
has a wave front set which is non- vanishing only in causal directions. Hence, thanks to the above 
cited theorem, we can conclude that the wavefront set of is empty or, equivalently, that it is 
smooth. □ 

Remark. A more detailed analysis of the scalar case and of the preceding theorem in particular 
can be found in theorem III. 2. 2. 7 of [HalO]. 



4.2 The Dirac Field Case 

The case of a Dirac field can be studied following the general ideas of the previous subsection, 
that is, we first construct a suitable state for the boundary theory and then prove that the 
natural bulk counterpart it induces is of Hadamard form. 

Proposition 4.2. The single-spinor two-point functions defined as 

toJ(*l,*J)= J 1 + \^p k *l(k, 9, y>)^(fc, 0, <p) dk rfS 2 

RxS 2 

are distributions on (5(Q~) (g) Via the double-spinor two-point function 

M*l,l © *1,2. *2,1 *S )2 ) = tO^(*i,i, *J >2 ) + W+(*2,l, ^ f i, 2 ) 

they unambiguously determine a quasifree state hp : — > C which is invariant under the 
automorphic action of the full boundary symmetry group Gq- . Furthermore, Ap fulfils the KMS 
condition at inverse temperature (3 with respect to the v— translations on 
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Proof. The two-point functions ro^ are linear and positive functional and, hence, they induce a 
unique quasifree algebraic state by well-known results (see chapter 5 in [BrRo96v2]). To verify 
the KMS condition and since ro 1 * 1 are manifestly invariant under ^-translations, we can basically 
repeat the proof for the scalar case: Let us set = VP^i ® ^| 2 > ^ ^ {!> 2} and let us introduce 

3(t) = tt)^(*i, a t (* 2 )) , &{t) = tt)^(a t (*i), * 2 ) , 

where at(^i)(v, 9, tp) = ^i(v — t, 9, tp). Both #(£) and <8(t) are bounded and continuous functions 
which can be seen as the Fourier-Plancherel transforms of two square-integrable functions $(E) 
and $(E) that satisfy the condition 

${E) = e^ E &(E) , 

which is a restatement of the KMS condition. Finally, invariance under the full automorphic 
action of Gq- follows as in the scalar case from the results of [DMP07], since ro^ are defined 
via invariant integrals. □ 



Remark. Although we do not prove it explicitly, we presume that one can show the just found 
KMS states on 25(Q _ ) to be unique and weakly converging to the unique (and pure) ground 
state with respect to ^-translations in the limit j3 — > oo. 

As in the scalar case, we can use proposition 3.4 to define an algebraic state A^ : 25(M) — > C 
induced by the boundary counterpart. Yet, before analysing the outcome, we have to introduce a 
further feature of a Dirac field theory, namely, the concept of charge conjugation: this operation 
is implemented by the *-automorphism a c : 23(M) — > B(M) defined as a c (B(f)) = B(f c ) for 
all / = f l e / 2 . Here, f c = C' l J 1 ® J 2 C where C G 5L(4,C) is the matrix which fulfils 
CC = I4 and 7 a = C^/ a C~ l for all a G {0,1,2,3}. Accordingly, it turns out that a quasifree 
state A : 'B(M) — > C is charge conjugation-invariant, i.e. Aoa c = A, if the associated two-point 
functions ro ± fulfil m ± (g^ c ,p c ) = to^(f,g) for all / € D(DM) and for all g G V(D*M). With 
this in mind, we proceed to prove the main result of this section. 

Theorem 4.2. The state kf : 25 (M) -> C defined as 

Af = Ap o %D . (29) 

possesses the following properties. 

a) The state A^ 1 is quasifree and charge conjugation-invariant. 

b) The two-point function tOM,/3 of A^ is an element of V (D M M D* M) of Hadamard form. 

Proof. Let us focus on a). As A^ is already quasifree, A^ is manifestly quasifree as well. The 
charge conjugation- invariance can be checked by explicitly computing that the single-spinor 
two-point function satisfies the necessary identities, see section III. 4 in [HalO] for details. 
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Let us now focus on b) and show that ro^- ^ G D'(DM M D*M) to start with. This result 
descends from the continuity of the causal propagator S with respect to the appropriate topolo- 
gies and from the continuity of the Fourier transform with respect to the relevant L 2 -norms. 
Since the modes associated to a solution of the Dirac equation display the same regularity as 
the modes associated to the Klein-Gordon equation conformally coupled to scalar curvature, we 
can proceed as in the proof of theorem 4.1 to obtain 

q(n) 

\\k n $ f \\ L2 <J2\\D j f\\L~, 

3=1 

where ^/ = T£>(Sf) and the notation is the same as the one in theorem 4.1. This part 

of the proof can then be finalised noticing that ro^^^, ^/ 2 ) are controlled by the L 2 norms 

||*/J|l»= / ** fi (h,6,4>)y fi (h,6,<l>)dvd§ 2 , i€{l,2}. 

RxS 2 

We are left with the verification of the Hadamard property. To this end, thanks to the results 
in [Hol99, KraOO], it is sufficient to ensure that the wave front sets of the distributions tu^- ^ fulfil 
the requirements of definition IV. 1 in [Hol99]. Furthermore, the charge conjugation-invariance 
of A^f implies that ro^ ^(/, h) c = tt)^ p(h}°, f^ c ). Hence, it suffices to analyse the wavefront set 
of ro~ and we have to prove that 

WF(xv M ^) = {(*, y, k x , k y ) G T*M 2 \ {0} | (x, k x ) ~ (y, -ky), k x > 0} . 

As ro^ f/3 = ro^ o (TpS (8> r^S 1 ), where S is the causal propagator and the projection of 
solutions to the horizon, this equality can be proven in the same way as the counterpart for 
the scalar field presented in [DMP09, PinlO]. The wave front set of this composition can be 
estimated along the lines given in the proof of theorem 4.2 in [Mo06] where theorem 8.2.13 in 
[Ho90] is shown to be applicable. □ 



4.3 The Bulk Remnant of the Boundary KMS Condition 

Our next goal is to discuss further physical properties of the states oj 1 ^ and Ajf , particularly, 
we would like to analyse the bulk remnant of the exact KMS condition on the boundary. As 
a starting point, we write the two-point functions of the states introduced in the preceding 
subsections in a more manageable form using the bulk mode-decompositions of the scalar and 
the Dirac field quantities. In the former case, the specific form of the causal propagators written 
at the end of section 2.2 entails that the two-point function of ui^f can be written as 

2A x ' X ' y ' V) ~ (2^a{r x )a{r y ) J \ 1- er^ + - 1 J & * * ' 
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where T\(r) are the modes discussed in and after proposition 2.1. In the Dirac case, we can 
make use of the analysis in section 2.3 to write the single-spinor two-point function defining 

as 

rv + M A^,x,T y ,y)= ^ e _ 0k + ^ -p— dk, 

R 3 «e{3,4} ie{i,2} 

where ^ l are the modes introduced in proposition 2.3 and appendix B. 

Notice that, although in the massless case a Bose-Einstein or a Fermi-Dirac factor at inverse 
temperature (3 are present in the integral kernel of the two-point functions, the bulk spacetimes 
are not invariant under time translations. Hence, the induced bulk states of the massless theories 
do not satisfy any KMS condition. However, a KMS-like condition can be found exploiting (2) 
and the conformal invariance of the considered scalar and Dirac theories. Particularly, we 
can relate each solution of the relevant equations of motion in a FRW spacetime with flat 
spatial sections to a counterpart in Minkowski spacetime by a rescaling with a suitable power 
of the scale factor a(r). Consequently, the induced states in the FRW spacetime can now be 
read as conformally rescaled genuine KMS states at inverse temperature (3 with respect to the 
Minkowskian time translations. This allows to recover a natural thermal interpretation. 

In the case of massive fields, the situation is more involved and a nice and simple analysis as 
the one above is not possible. Yet, it would be nice if the induced bulk states still enjoyed some 
properties which could be seen as the remnant of the exact KMS on the boundary. It seems 
natural to approach the issue by estimating how much the regularised two-point functions deviate 
from those of the conformally invariant case. We shall analyse only the case of a scalar field, and 
will briefly comment on the Dirac case afterwards. Hence, we are concerned with the difference 



(w/3,2 - ^oo,2)(r x ,x,T y ,y) 



1 f ( T k (r x )T k (r y ) + T k {T x )T k {T y ) \ 3 

(27rfa(T x )a(T y ) J V e^-l J 

»3 



where the dependence on the mass affects only the form of the modes T k . The fact that these 
are constructed as a convergent perturbation series which fulfils the initial condition of being 
asymptotically of the form s/2k e lkr , i.e., resembling massless positive frequency modes, shall 
allow us to compare massive and massless theories. For the sake of simplicity, we discuss such 
a relation only for a specific observable, the Wick-regularised squared scalar field :0 2 :. In the 
massless and, hence, conformally invariant case we obtain 

«-^)(0 2 (r„f)) 1 



In the massive case, in order to compute the expectation value of the same observable, it suffices 
to notice that the perturbative construction and initial conditions on T k imply 



(wwi'-s)!^^- 1 ) 
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This inequality can be used to control the deviation from thermal equilibrium, as it can be 
directly seen from the analysis of the expectation value of : 4> 2 : , viz. 

l«'(^-«J0(^)l^ + I ^(^-i). 

Although we have only discussed the influence of the mass on a rather special observable, the 
outcome of the same procedure in other cases such as, for example, the balanced derivatives 
introduced in [BuOjRo02] for the flat case and generalised to curved spacetimes in [SchVe08]: 

9(n,/3,x) = Um(Vni • ■ ■ ^n H (f>)(^>i3 ~ Woo)(r + £,£,r - £,x), 

can be computed in a similar way. In the massless case, the expectation values of these observ- 
ables turn out to be fixed and simple functions of the temperature. However, in the case of 
massive fields, the fact that the thermal nature of these states holds only in an approximated 
way manifests itself as a modification of the relations existing between different G - so-called 
transport equations - due to the appearance of central terms, i.e. sources for the transport 
equations, depending on both the mass, a(r), and its derivatives. 

A similar analysis as the one above is certainly possible in the case of Dirac fields. However, 
the observable which seems to be the natural one in this case, the Wick-square regularised 
with the suitable single-spinor two-point function of A^, is, in contrast to its scalar counterpart, 
not a good thermometer since its expectation value in the massless case vanishes. 

4.4 Relaxing the KMS Condition on the Boundary 

In the previous subsection we have found that, in the massive case, the induced states, which 
fulfil an exact KMS condition on have an approximate thermal interpretation in the bulk 
which is quantitatively controlled by am. Hence, one can interpret these states as states which 
are in perfect thermal equilibrium at a = 0. From this point of view, a natural generalisation 
of the states constructed above are states which are in thermal equilibrium "at an instant of 
time" where a = do > 0. To wit, we can define a new set of states oj^ ao and A^ aQ as the (charge 
conjugation-invariant) quasifree states induced by the two-point functions 



UJoR (t XT if) — - f ( Tk ^ Tk ^ _ + T k(Tx)T k (T y ) \ t k(£-y) ^ 



R3 

and 



*° + M,p,ao(Tx,X,Ty,y) = 2^ a 7m 2 ' 2 + o /ux 2 2 dk - 

These states are still Hadamard as one can prove by trivial modifications of the relevant proofs 
for the case qq = 0. However, the corresponding states on the boundary algebra certainly fail 
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to fulfil an exact KMS condition. Nevertheless, from a physical point of view, one can expect 
that the states with ao > are better suited to describe realistic thermal situations in the early 
universe than the states with ao = 0. 



4.5 Approximated Thermodynamics 

As last point, we briefly sketch how approximated thermodynamic relations hold for the above 
introduced states, both for ao = and for ao > 0. In order to support the interpretation of these 
states as approximated equilibrium ones, we analyse whether and to which extent the usual laws 
enjoyed by any theory with a sensible notion of thermodynamics hold for certain expectation 
values. For the sake of simplicity, we consider the scalar field states ^^ ao - Furthermore, in order 
to avoid potential complications due to well-known anomalies as the trace anomaly, we shall 
consider the difference and, hence, define 

(A) fl „=ufi ao {A)-ug{A). 

As we wish to formulate thermodynamic laws with respect to the evolution generated by 
e = dt where t is the cosmological time, we introduce, following [Dix78, BuOjRo02], the vector 
representing the effective inverse temperature as 



{<p 2 (x))f3,a ' 

Since both and are of Hadamard form, (4> 2 (x)) p <afj is a well-defined smooth function on 
M. Moreover, this definition allows us to decompose the expectation values of the stress-energy 
tensor T^ v evaluated once more with respect to w^ ao — as 

where Q represents the heat and P the pressure. According to [Dix78], we can additionally 
introduce the vector representing the thermodynamic entropy as 

S(p c y = Q(f3 c )f3 c e» . 

With these definitions for /3 C , Q, P and S at hand, we can discuss potential thermodynamic 
relations among them. In fact, one can show that the following laws hold near a = and/or for 
small m. 

0th law: f3 c is constant on the Cauchy surfaces at fixed cosmological time. 
1st law: (Tp^p^ao is covariantly conserved. 



2nd law: The entropy production is positive, that is 

4 
3 
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3rd law: The minimum of the entropy function is reached in the limit (3 C — > oo. 

The most difficult relation to prove is the second one. In the massless case, the entropy pro- 
duction is vanishing because (3 C turns out to be exactly a conformal Killing vector and (T)p jao 
vanishes. In order to address the massive case with ao = 0, we can make an expansion of the 
relevant quantities in powers of m and finally restrict our attention to the region near Q 5- . Per- 
forming such an approximation, we observe that /3™ computed for a mass m, differs from the 
massless one /3° via 

/3 c m -/3 c = (- Cl mV) + O(m 3 ), 

where c\ is a positive constant. Furthermore, the first order in m 2 of {T)p m is equal to — C2^- 
where C2 is another positive constant. Finally, we notice that \(T^ - \Tg^) PA0 \ < Ca~\ and 
that (T^ v — ^Tg^ u ) p iao ¥ (nficu) vanishes because j3® is a conformal Killing vector. With these 
observations in mind we can compute a 2 VS up to 0(m 3 ) terms and notice that it tends to a 
positive constant towards Hence since a 2 is positive and since a 2 VS is smooth in M, VS 
must be positive near 9 + . In the case ao > 0, we have only been able to verify the positivity of 
the entropy production by numerical methods, as the occurring integrals can not be evaluated 
analytically any more. 



5 Conclusions 

We have studied the classical and the quantum behaviour both of a real, massive, conformally 
coupled scalar field and of a massive free Dirac field in a large class of Friedmann-Robertson- 
Walker spacetimes {M^gpRw) with flat spatial sections. These spacetimes are singled out by 
criteria which have both a physical and a mathematical origin. Particularly, we have constrained 
the dependence of the scale factor a(t) on the cosmological time t to be such that the range 
of the conformal time r is unbounded on the negative real axis. Moreover, we have demanded 
that the dependence of a(r) on r is of the form (4), which entails that the corresponding 
Universe underwent a phase of either exponential or power-law expansion in its early stages. 
From a geometrical point of view, our requirement on a(t) entails that the underlying manifold 
can be endowed with a boundary Q~ which is either a cosmological horizon or a Big Bang 
hypersurface, and mathematically a null differentiable manifold of codimension 1. On 9~ it 
is possible to intrinsically construct both a scalar and a Dirac field theory whose associated 
algebras of quantum observables are sufficiently big to contain the bulk counterpart via an 
injective *-homomorphism. The motivation to introduce this injective homomorphism lies in 
the particular topological structure of which allows to define an exact notion of a quasifree 
KMS state up for the boundary theory with respect to the rigid translations along the complete 
null geodesies generating Furthermore, the above-mentioned injection allows to pull-back 
up to the bulk, yielding a quasifree state which, both in the scalar and in the Dirac case, has 
various physically and mathematically relevant properties. To wit, such a state always fulfils 
the Hadamard condition and, when the mass vanishes, also the KMS one. 
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Prom a physical point of view, our results have two immediate possible applications: on the 
one hand, we can now extend the analysis of [DFP08] and look for stable, homogeneous, and 
isotropic solutions of the semiclassical Einstein equations driven by quantum fields of higher 
spin. This has been done in [HalO, DHMP10], and interesting results have been found. On 
the other hand, the states we have introduced are both Hadamard and possess an approximate 
thermodynamic interpretation also for massive fields. Hence, one can expect that they provide 
a natural starting point to analyse the cosmic neutrino background on a genuine quantum field 
theoretic level. It can be expected that the deviations from thermal equilibrium give rise to 
modifications of the standard picture, in contrast to the cosmic microwave background, which 
relies on a well-defined thermodynamic interpretation, as photons are massless. We hope to 
return to this issue in the near future. 
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A On the Definition of Dirac Fields 

The role of this section is to provide a short introduction to the mathematical structures which 
are needed to define the notion of Dirac fields on a four-dimensional globally hyperbolic space- 
time. Our goal is not to provide an exhaustive analysis but rather to offer to a potential reader 
the chance to understand the content of this paper without necessarily resorting to other long in- 
troductory manuscripts. Hence, this appendix shall summarise what has already been presented 
in [Sa08, DHP09, HalO]. 

In this paper we only consider Friedmann-Robertson- Walkers spacetimes M, which are 
four-dimensional, Hausdorff, connected, smooth, simply connected, oriented and time-oriented 
globally hyperbolic manifolds endowed with a Lorentzian metric, whose signature is chosen as 
(— , +, +, +). Therefore, we can introduce two important ingredients: 

• The spin group Spin(l,3) as the double cover of 50(1,3), i.e., there exists the following 
short exact sequence of Lie group homomorphisms: 

{e} — ► Z 2 — ► Spin(l, 3) — ► 50(1, 3) — ► {e} , 

where {e} stands for the trivial group, whereas Z 2 = {±1} is the cyclic group of order 2. 
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• The frame bundle associated to the tangent bundle TM, that is, the principal bundle 

FM = FM[SO (3,l),7r / ,M], 

where S0q(3, 1) is the component of the Lorentz group connected to the identity. Notice 
that FM is build from the disjoint union \_\ x F X M, where F X M is identified with the typical 
fibre SOo(3, 1) and where ir' : FM — > M is the projection map. 

We are now in the position to introduce the main geometric structure at the heart of the 
construction and of the analysis of Dirac (co)spinor fields. 

Definition A.l. Given an oriented and time oriented spacetime M , a spin structure is the 

pair (SM,p) where SM = SM[Spino(3, 1), 7r, M] is a principal fibre bundle over M with the 
identity component of the spin group as typical fibre. Moreover, p is a smooth equivariant bundle 
morphism from SM to FM, that is, the following two conditions hold: 

1. p is base point preserving, such that 

ir' O p = 7T, 

2. p must be equivariant, i.e., calling R^ and R\ the natural right actions of Spino(3,l) on 
SM and of SOo(3, 1) on FM respectively, we require that 

poR x = R A op, VAG5O (3,l), 

where A = 11(A), being IT is the surjective covering from Spin(3, 1) to SO(3, 1). 

For the scenario we are interested in, the fact that M is four-dimensional and globally hyperbolic 
not only assures that a spin structure exists, but this structure is also unique up to equivalence, 
since M is a simply connected spacetime [Ge68, Ge70] (c./., [Dim82] for the definition of equiva- 
lence in this context). The kinematically allowed configurations can now be readily defined out 
of the objects at hand. 

Definition A. 2. We call Dirac bundle of a four dimensional Lorentzian spacetime M the C 4 - 
bundle DM = SMxyC 4 associated to SM with respect to the representation T = ,0 )©£)( '2) 
of SX(2,C) ~ Spino(3, 1). This is the set of equivalence classes [(p,z)], where p £ SM, z G C 4 
and equivalence is defined out of the relation 

(pi,zi) ~ (P2,z 2 ), 

if and only if there exists an element A of SL(2, C) such that Ra(pi) = P2 and T{A~ 1 )z\ = z-i- 
The global structure of DM is that of a fibre bundle over M with typical fibre C 4 , and the 
projection map ttd is traded from the one of SM, namely, V [(p, z)\ 6 DM, it holds 

kd[(p,z)] = 7r(p). 
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Furthermore, if we endow C 4 with the standard non- degenerate internal product, we can construct 
the dual Dirac bundle D*M as the C A *-bundle associated to SM requiring that {p\,z\) and 
(Ra(pi), z*T(A)) are equivalent, where * denotes the adjoint with respect to the inner product 
on C 4 and elements of C 4 * are understood as row vectors. Consequently, the dual pairing of C 4 
and C 4 * extends in a well-defined way to a fibrewise dual pairing of DM and D*M. 

According to this definition, one can introduce the following notions. 

• A Dirac spinor is a smooth global section of the Dirac bundle, i.e., if) € E-(DM). Since 
DM is trivial due to M being four-dimensional and simply connected, is (diffeomorphic 
to) a vector- valued function ip : M — > C 4 . 

• We call Dirac cospinor a smooth global section of the dual Dirac bundle, namely, ip' G 
£(D*M). On a FRW spacetime, V' : M -> C 4 *. 



B Perturbative Construction and Analysis of the Classical So- 
lutions 

In order to have full control on the behaviour of the classical solutions of both the scalar and of 
the Dirac Cauchy problems, we must discuss the particular form of the functions Tj, and \k which 
are present in the expansions (7) and (16), respectively. Since the analysis of the scalar case 
has been already performed in [DMP07, DMP09, PinlO], the aim of this section will be to only 
address the Dirac case. If we focus on (16), then the construction of each \k can be performed 
along the same lines pursued in the scalar case. Moreover, we only discuss the asymptotic de 
Sitter case specified by 5 = in (4) and briefly mention the modifications necessary to treat the 
asymptotic power-law case, i.e. 5 > 0, afterwards. 

In other words, we first consider (16) in the de-Sitter Universe and, then, we solve the general 
case by means of a convergent perturbation series. Hence the starting point is (16) written as 

^ + V (k,T) + V(T)^Xk(r) = 0, (30) 

where V (k, r) = k 2 + - ij^j 4? whereas V(t) = k 2 + a(r) 2 m 2 - m'(r)m - V (k, r). Notice, 

that, according to (4) with 7 = —H^ 1 , the function V(t) is either or at least of order 0(r~ 3 ). 

Since we would like to treat V(t) as a perturbation potential, we obviously start from the 
unperturbed solution. If V(t) = we are in a pure de Sitter spacetime and two solutions which 
satisfy the desired initial conditions are 



XL = A^H^i-kr) X 2 2 = B^H^i-kr), (31) 



Xk,2 

where J#> and are the first and second Hankel functions with 
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Furthermore, in order to satisfy the desired normalisation (17) we have 



A = =— , B = 



2k V2k 

Since we have a full control of the solutions of (16) in the cosmological de Sitter spacetime, we 
can now revert back our attention to (30) and we shall simply proceed along the same lines of 
[DMP07, DMP09] for the real scalar field. Hence we write a formal solution as a Dyson-Duhamel 
series: 

X^o) = x£,i(r ) + ^(-ir / dn / dr 2 ... / dT n Y[S k (n-i,ri)V(T^Xk,i(r n ), (33) 

ra=] ?— 1 

a 1 -oo -oo -oo % — 1 

where tq = t and where Sk{r,T r ) is the retarded fundamental solution of (16) in cosmological 
de Sitter spacetime: 



S k (r,r') = -^( X m (t)x° i2 (t') - X^iW^r)). 



Notice that, being antisymmetric, Sfc(r, r) = 0, while the time derivative is a conserved quantity 
of the underlying dynamical system and thus we can simply set dSk{ d l y) = 1. Of course 

T=T ! 

(33) would be rather useless if we were not able to prove its convergence and, to this avail, we 
need to provide suitable estimates for both x k j an d for the involved integrand. 

Let us start from Sk(r, r') and, since we would like to have an estimate for Xfc( T o) which 
guaranties square integr ability in k and at most polynomial growth for large k, we shall give two 
estimates, for small and large k. Let us start analysing the large k behaviour. We observe that, 

being Hi 2 \z) = Hg\z), we can make use of the integral representation present in formula 9 in 
8.421 of [GR95], in order to express both x^j m a m ore manageable form. We notice that the 
integral representation can be uniformly estimated in order to give 

\xh\ ^ c »l> ( 34 ) 

where C v is some positive constant. With this result we obtain, by direct inspection that 

\S k (T,T')\<C± (35) 

holds uniformly in r and r' for some positive constant C. The latter will be suitable for analysing 
the convergence of the perturbative series for large value of k. We proceed now to derive another 
estimate in order to control the behaviour of |xfc| for small values of k. We notice that S k (r, r') 
can be rewritten in terms of the Bessel J„ functions in the place of the Hankel one, to obtain 
that 

\S k (r,r')\ < ClV^?\j u (-kr)J^(-kr') - J v {-k T ')J-u(-kT)\ , 
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where C v is again a certain positive constant. We can now make use of the recursive relations 
to rewrite J- V {z) in terms of the J- u +i and its first derivative (formula 1 in 8.472 in [GR95]). 
Now both Bessel functions can be estimated analysing their integral representation present in 
formula 5 in 8.411 in [GR95]. In that way we obtain the following estimate valid for |r| < |r'| 

\S k (r,r')\<C(\T'\ + k\rf), (36) 

where C is some positive constant. 

With these estimates in mind, we can now analyse the convergence of the series (33). From 
(36) we obtain that the following uniform estimate holds 

\Xk(r) - xli(r)\ < f exp / C (|n| + k |n| 2 ) - 1 ] \V{n)\ dn sup | X ;1 (r 2 )| , (37) 

V OO / T2e(-0O,T) 

while, from (35), the series (33) can be uniformly estimated by 

\Xk(r) - x2,i(r)| < f exp / £ - 1 ) \V(n)\ dn sup | X ° M (r 2 )| . (38) 

V -oo / T 2 e(-oo,r) 

If |V(t)| < C|r| 3+e , both integrals present in the preceding inequalities are finite. 

We proceed now to briefly discuss the approximation we shall use for the time derivative of 
Xfc(r). Since the employed procedure is the same as the one written above, we simply summarise 
the result. We obtain that also the series for the first order derivative of Xk in the time variable 
is uniformly bounded by a convergent series, namely 



T 



< 



<C I W 2 {k,T,T 2 )\V{T 2 )\dT 2 exp f Wi(n,fc) |y(ri)|dn sup |x°,i( r 3)l • (39) 

J J T3£(-oo,t) 

— oo — oo 

where Wi(T\,k) is a positive real function greater than \S)~(t, ti)| uniformly in r. Hence, ac- 
cording both to (36) and to (35), W\(ti, k) can be either C ( | r 1 1 + &|ti| 2 ) or k" 1 . Furthermore, 
W 2 (k,T,T 2 ) is another positive real function that controls uniformly the derivative of the propa- 
gator, that is \-§pSk(r, t 2 )\ < W 2 (k, r, t 2 ). Proceeding as above we observe that W 2 can be either 
C(l + l^|), which is an appropriate bound for large values of k, or C(l + + \kr 2 \ + \kr 2 \ 2 ), 
which is better suited for small k. We summarise the discussion as: 

Lemma B.l. For sufficiently large values of \t\, the series (33) for both Xk an d its first time 
derivative d T Xk is uniformly convergent i/|V(r)| < C|r| 3+e . Furthermore, for fixed t, both Xk( T ) 
and d T Xk(T~) are contained in L 2 ([0, k±], k 2 dk) for any k\ > 0. For large values of k both Xfc( r ) 
and its first time derivative are at most growing polynomially. 
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Proof. The convergence of the series descends easily from the discussion presented above. The 
local square integrability shown by Xk descends from (37) while the polynomial boundedness for 
large k from (38). The very same properties shown by d T Xk are obtained from (B). □ 

As anticipated, we now briefly mention how analogous results can be obtained in the asymp- 
totic power-law case. In this case, the solutions can be constructed as a perturbative series 
around the solutions of the massless, conformally coupled Klein-Gordon equation. Hence, 
Vq{t) = 0, whereas V[t) = ia'{T)m + a 2 (r)m 2 . Suitable estimates for V(r) follow trivially 
from (4), consequently, the convergence and fc-regularity properties of the perturbation series 
can be discussed as above. We refer the reader interested in further details to the proof of lemma 
II.4.1.5 in [HalO]. 
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